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Abstract 

We make a brief review of the Kramers escape rate theory for the 
probabilistic motion of a particle in a potential well U(x), and under the 
influence of classical fluctuation forces. The Kramers theory is extended 
in order to take into account the action of the thermal and zero-point ran- 
dom electromagnetic fields on a charged particle. The result is physically 
relevant because we get a non null escape rate over the potential barrier 
at low temperatures (T — » 0). It is found that, even if the mean energy 
is much smaller than the barrier height, the classical particle can escape 
from the potential well due to the action of the zero-point fluctuating 
fields. These stochastic effects can be used to give a classical interpreta- 
tion to some quantum tunneling phenomena. Relevant experimental data 
are used to illustrate the theoretical results. 
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1 Introduction 

One of the most useful contributions to our understanding of the stochastic 
processes theory is the study of escape rates over a potential barrier. The 
theoretical approach, first proposed by Kramers [Q, has many applications in 
chemistry kinetics, diffusion in solids, nucleation [2], and other phenomena [3]. 
The essential structure of the escape process is that the bounded particle is 
under the action of three types of forces: a deterministic nonlinear force with 
at least one metastable region, a fluctuating force whose action is capable of 
pushing the particle out of the metastable region, and a dissipative force which 
inevitably accompany the fluctuations. 
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Figure 1: Metastable potential with a barrier height AU = U(xb) — U(x a ), a 
local minimum x a and a local maximum x\> (top of the barrier) . 



In this work we describe the escape rates of a particular model: a classical 
charged particle moving in the metastable potential shown in the Figure 1, and 
under the influence of the fluctuating electromagnetic radiation forces commonly 
used in Stochastic Electrodynamics (SED)0]|H]. The fluctuating fields postu- 
lated in SED are classical random fields, with zero mean but nonzero higher 
moments. The spectral distribution of this radiation can be expressed as a sum 
of two terms 
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The first term is the zero-point radiation contribution to the spectral distribu- 
tion. It is independent of the temperature and is Lorentz invariant. The second 
term in Q is the blackbody radiation spectral distribution, responsible for the 
temperature effects on the system. 



2 Properties of the harmonic oscillator motion 
under the action of a thermal and zero-point 
radiation 

The zero-point radiation (first term in has a mean energy fno/2 associated 
with each mode of the electromagnetic fields, and is responsible for the most 
important features of SED. With this zero-point radiation postulated, several 
phenomena associated with the quantum behavior of the microscopic world can 
be explained on classical grounds. Many examples can be found in the reviews 
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The potential U(x) will be approximated by a harmonic oscillator in the 
region of the potential well [x « x a ) so that (see Figure 1) 

U(x)~U(x a ) + -muj 2 a (x-x a ) 2 , (2) 

where u a is the natural frequency of the oscillator. 

The dynamical behavior of a harmonically bounded charged particle has been 
extensively studied in the context of classical SED. It is found that the zero-point 
radiation maintains the stability of this system. We shall use the statistical 
properties of the harmonic oscillator in order to understand, classically, the 
escape rate at very low temperatures. We give below a brief review of the 
harmonic motion under the action of the random electric fields characteristic of 
SED. 

The nonrelativistic motion of the charged particle (charge e and mass m) 
near the bottom of the potential well (see Figure 1) is governed by the equation 

2e 2 - 

r< = -m^ +eE x (t), (3) 

where £ — x — x a , and E x (t) is the x component of the random electric field. 
The term proportional to £ is the radiation reaction force. The electric field is 
such that (E x (t)) = and 

{E x {t)E x {Q)) = — / dwp(w,T)cos(wf), (4) 
Jo 

where the spectral distribution p(u), T) was introduced in the equation ft). The 
radiation reaction force can be approximated by [5] 

2e 2 - 

_£~-m 7 £, (5) 

where 7 = 2e 2 w 2 /3mc 3 . Moreover it is verified that 7«w a . According to these 
approximations one can show that the average energy (e) of the oscillating charge 
is such that 

<e> = + \m£(e) = ^ coth (|^) ee D(T), (6) 

where we have introduced the function D(T) in order to simplify our nota- 
tion. Notice that the average energy depends on the temperature and on the 
oscillatory frequency ui a . 

The result © is well known j^j. The average energy (e) becomes equal to 
ksT in the high temperature limit (kgT 3> huj a ), and is non zero when T = 0. 
Actually D(T) -> huj a /2 as T -> 0. Notice that D(T) depends on h. We can 
show that the Planck constant comes from the intensity of the zero-point field 
E x (t) that appears in We recall that hui a /2 is the value of the ground state 
energy of the harmonic oscillator in quantum mechanics. This result, obtained 
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within the realm of SED, differs from the usual null result of ordinary classical 
physics because the zero-point fluctuations are taken into account. 

It is quite natural to use the average energy © in the calculation of the 
Kramers escape rate of a potential well. We shall see that the consequence 
of the new form of the average energy is a non-vanishing escape rate even if 
T — > 0. Most authors do not mention the classical zero-point fluctuations and 
use the quantum mechanical formalism to interpret the non null escape rate 
as a tunneling through the classical forbidden region of the barrier. We shall 
see that the zero-point fluctuations allow the escape over the potential barrier 
even if the mean energy of the particle inside the barrier is much less than 
AC/ = U(xb) — U(x a ). In the classical mechanics context the escape would be 
impossible without the action of the zero-point fluctuations. For simplicity we 
shall take U(x a ) — in what follows. 



3 The escape rate over the potential well 

We shall use the approach of Chandrasckhar ^Uj , based on the Kramers theory. 
The physical system considered by Chandrasekhar is a particle moving under 
the influence of a fluctuating force, and a potential U{x) that has a metastable 
region (see Figure 1). The motion of the particle is governed by a Langevin 
type equation 

mx = -iwyx -U'(x) + F(t), (7) 

where —rwyx is the dissipative force and F(t) is the fluctuating force which is 
characterized by the average (F(t)) — 0. The average energy of the particle 
within the potential well, that is x < Xi, is assumed to be given by 

(jX 2 + U(x)} = k B T, (8) 

in the high temperature limit. 

It is possible to show that the Langevin equation J7J) leads to a phase space 
Fokker-Planck equation given by JU| 

-^r + -^--U'(x)—= 1 W + yp—+m 1 k B T—-, (9) 
at m ox op op op z 

where W = W(x,p,t) is the probability distribution in phase space. Notice that 
the left hand side of the above expression is equivalent to the Liouville equation. 
The right hand side appears as a consequence of the fluctuating and dissipation 
forces. For low temperatures, the Fokker-Planck equation @ is not valid. As 
we have mentioned in previous section (see the equation ©), the factor fcgT in 

the last term of © must be replaced by D(T) = ^ coth f at^r) ■ Therefore, 
we shall consider the following equation 



dW p dW ,dW „, dW ^^d 2 W 

— + -—- - U'(x)— = 7 W + 1P — + m 1 D{T)— r 
at m ox op Op op 
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We shall sec that the equation I|1U|) will allow us to give an accurate description 
of the escape rate at low temperatures. 

In the Kramers theory, two quantities are essential to calculate the escape 
rate. One is the probability P(t) of finding the particle inside the potential well. 
This probability can be obtained from the phase space distribution, namely 

/oo rx b 
dp / dxW(x,p,t). (11) 
-oo J — oo 

The other important quantity is the diffusion current, j(xb), across the top of 
the potential barrier. The diffusion current in an arbitrary position x is defined 

by 

j(x,t)= f dp^W{x lP ,t). (12) 
Using (fTOf) , (fTT|> and (JI35J one can show that 

(> ' ,{l) - dp^-W(x b ,p,t) = -j(x b ,t). (13) 



dt J_ 00 m 

The escape rate k, regarded as the decay factor of the probability P(t), can 
be defined by the equation 

^ = -«P«). (14) 

The solution of the above equation is 

P(t) = P e- Kt , (15) 

where Pq is a constant that will be calculated later. On the other hand, con- 
sistently with the equations (fT^jl and lfH|) . one can define the escape rate as 



_ j{Xb,t) 



(16) 



P(t) ■ 

Therefore, according to the above theory we have 

W(x, Pl t) = Q(x,p)e- Kt , (17) 

where Q(x,p) satisfies the equation 

p 8Q ,dQ dQ ™d 2 Q 

— U (x) jQ - ~fp- mjD(T) — - = 0. (18) 

m ox op Op op z 



A physically interesting case could be 

p 2 /2m + U(x) 



Q(x,p) oc exp 



D(T) 



(19) 



however, this standard distribution leads to a situation in which there is no 
diffusion across the potential barrier at Xb- 
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Under the conditions of our problem, the equilibrium distribution can- 
not be valid for all values of x. Hence we shall consider a solution of the equation 
1181) in the following form 



Q{x,p) = CF(x,p) exp 



D(T) 



(20) 



where F(x,p) is an unknown function that will be determined below, and C is a 
normalization constant. Notice that the above expression is valid for the phase 
space motion near the bottom of the potential well (see Figure 1). The new 
form for Q(x,p), introduced in (|2()(l . requires a boundary condition on F(x,p) 
for x ~ x a , namely F(x,p) ~ 1. An alternative form for this boundary condition 
is PU 

F(x,p) -» 1, for x < x b . (21) 

Another physical hypothesis is necessary. The probabilistic motion near the 
top of the barrier (x ~ Xb) is also governed by a phase space distribution similar 
to (J20J), namely 



Q{x,p) = CF(x,p)exp 



^ + U(x b ) - \mu>l{x - x b ) 2 



because 



U(x) ~ U(xb) - -muil(x - x b y 



(22) 



(23) 



The use of D{T) in both formulas Q20f> and (|22|l is justified because we are 
assuming that the particle stays a long time in the potential well (see section 
2), and crosses the top of the barrier very quickly (uj b 3> oj a ). 

Since only a few particles can escape over the potential barrier, another 
boundary condition must be imposed on F(x,p), that is 



F(x,p) — > 0, for x ^> x b - 



(24) 



Notice that the boundary conditions l|21|l and l|24|) are simple hypothesis that 
can be justified on physical grounds. 

Substituting (1221) into l|18|) . we obtain the following differential equation for 
F(x,p), 



pdF 2 OF OF nfT s9 2 F 

m ox op op op* 



(25) 



Following Chandrasekhar we assume that F(x,p) = F (p — am(x — x b )) = 
F(y), where a will be obtained below. With the introduction of the variable y, 
we obtain the more simple differential equation 



, dF , d 2 F 

. (a _ 7)y _ =mTD(T) _ 



(26) 
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provided that the constant a is such that 



a — 7 

This equation for the constant a has the solutions 



7 / 7 



The single variable differential equation l|26|) can be integrated giving 



F = Fq exp 



(a - 7)y 



/2 



2m-yD(T) 



dy', 



(27) 



(28) 



(29) 



where Fq is a constant. One can see that only the positive root in (|28|l leads to 
a — 7 positive, so that F(y) naturally obeys the boundary conditions l|21|l and 
Q24|). Therefore, one can show that 



F(x,p) 



a — 7 



2TTm-/D(T) 



exp 



(a - 7)2/ 



'21 



2mjD(T) 



dy 1 . 



(30) 



Combining l|22[) and (|30|l we get for Q(x,p) the result 

Q(x>P) = C X [Z llZrnS X ( 31 ) 

(a-7)y' 2 " 



x exp 



+ Z7 (x fc ) - \mu}l(x - x h f 



D(T) 



exp 



2m-yD(T) 



dy 1 . 



The equation (13 1|> is valid only in the neighborhood of Xb- Inside the potential 
well (x ~ x a ) the approximate solution is (see section 2) 



Q{x,p) = C exp 



hi + \ muj \{x - x a f 
D(T) 



(32) 



Using the expression l|32|l . and considering the equations 115|1 and (|17() . 
we obtain for the constant P 



OO />oo 



— OO •/ —OO 



dp J dxC exp 
27rD(T) 



^ + \mu 2 a {x~x a y 
D{T) 



(33) 



The diffusion current across the top of the barrier is (see (|12|l and H17|) ~) 



dp—Q(x = x b ,p), 
m 



(34) 
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where Q(x = Xf,,p) is given by our expression (|31[l . From (|34[1 we get 



j(xb) = C 



a — 7 



27tto7D(T) 

dp — exp 

m 



exp 



C/(z b ) 



£>(T) 



(35) 



2mD(T) 



It is straightforward to show that 



j(x b ) = cJ^—-^D(T)exp 



The escape rate k, defined in (Tffi , becomes 
i(^fc) 



2tt 



a — 7 



■ exp 



<£?/ exp 



'D(T) 



(a - 7)y' 



2m^D{T) 



D(T) 



(36) 



(37) 






' AU ' 


J exp 


[ D(T)\ 



independently of the normalization constant C. Notice that the exponential 
factor e _Kt , present in j(xi,,t) and P(t), cancels leading to the result (|37|l . 
From the expression l|28|) for the positive root, it is possible to show that 



(38) 



Notice that in the low friction limit 7 — > 0, this expression is the simple for- 
mula K ~ (w a /27r) exp (— AU/D). It is very important to remark that, in this 
equation, the escape rate depends on the potential height AU, and on the 
parameters characterizing the particle motion inside the barrier, namely, the 

frequency uj a and the average energy (e) = D(T) = ^p- coth | 2 ^° T j . We recall 

that (e) = ksT when the temperature is high enough. 



4 Comparison with experimental data and con- 
clusion 

In order to illustrate, in a quantitative manner, the great analogy between the 
quantum tunneling description and our classical stochastic escape rate calcula- 
tion, the experimental results of Alberding et al. will be used. Notice that 



7 7 e 2 hui, 

— < — = 

ujb foe fnc 

so that the expression (|38|l can be written in the form 



«7 L = 1-^«1> (39) 
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Figure 2: The solid line gives the high temperature behavior of our theoretical 
prediction (formula 140(1 ). The intermediate temperature experimental data 
for the escape rate associated to the CO migration to a separated j3 chain of 
hemoglobin are plotted (Alberding et a!., 1976). The open circles corresponds 
to the escape rate without the zero-point fluctuations, namely the Arrhcnius 
formula k = ^ exp{~AU/k B T). 



According to Alberding et al., the beta-chain of hemoglobin ((3Hb) is bounded 
to the carbon monoxide CO from which it can be separated with a LASER. The 
rate of recombination can be obtained experimentally. The fraction N(t) of the 
molecules that have not been recombined with CO is measured as a function of 
time. Then, the time r, necessary to reduce N(t) to 75% of its original value, 
is determined. It is assumed that this recombination is a passage through the 
potential barrier and a good estimate of the escape rate is k = 1/r. 

This experimental procedure can be repeated for different temperatures T. 
The result for k(T) is indicated by the experimental points (black dots) in the 
Figures 2 and 3, obtained by N.R. Alberding and collaborators. We shall see 
that these experimental data are very well described by the formula (|40[1 . 

We have adjusted the values of cu a and At/ so that the experimental data 
and the formula (|40|l arc in good agreement. The values obtained are 

^ - 2.53 • KT 3 eV, 
2 

AU = 6.68 • 10 _2 eV. (41) 
Notice the impressive agreement between the classical theory with zero-point 
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Figure 3: Experimental data for the escape rate associated to the CO migration 
to a separated (3 chain of hemoglobin (Alberding et al., 1976). The solid line 
is our theoretical result (see formula l|4U|l ). The open circles correspond to the 
escape rate without the zero-point fluctuations (Arrhcnius formula). 

radiation and the experimental data. This is more clearly seen in the Figure 
3. We conclude that the particle can escape from the potential well at T — > 0, 
despite the fact that the barrier height AU is much bigger than the particle 
mean energy fiw a /2 inside the well (2AU/hui a « 26). 

It is interesting to recall that Alberding et al. have obtained AU using a 
conventional quantum mechanics calculation. They have found a value for AU 
which is in semi-quantitative agreement with the value obtained by us. However, 
the frequency uo a was not obtained by Alberding et al.. We want to stress that 
the frequency uj a gives relevant information about the potential well (see (|2|). 

The classical stochastic interpretation of the zero temperature escape rate 
is that the zero-point fluctuations provide enough energy so that the particle 
can go over the potential barrier. Since the particle is subjected to both the 
fluctuation and the dissipation processes associated with the radiation bath, the 
energy is not a constant of the motion. Therefore, particles that are initially 
inside the potential well can escape and be detected at points x c 3> (see 
Figure 1), with a fluctuating energy e(t) < AU, contrary to the criticism of 
Baublitz concerning the SED type calculation • Therefore, the classical 
escape rate calculation presented in our paper gives results entirely analogous 
to the quantum tunneling description, provided that the electromagnetic zero- 
point fluctuations are included in the calculations. 
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